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ABSTRACT 

Recent observations have shown that the direct-impact Algol systems U CrB and 
RS Vul possess gas located outside of the orbital plane, including a tilted accretion 
disc in U CrB. Observations of circumstellar gas surrounding the mass donor in RS 
Vul suggest magnetic effects could be responsible for deflecting the accretion stream 
out of the orbital plane, resulting in a tilted disc. To determine whether a tilted 
disc is possible due to a deflected stream, we use three-dimensional hydrodynamic 
simulations of the mass transfer process in RS Vul. By deflecting the stream 45° out 
of the orbital plane and boosting the magnitude of the stream's velocity to Mach 30, 
we mimic the effects of magnetic activity near LI. We find that the modified stream 
parameters change the direct-impact nature of the system. The stream misses the 
surface of the star, and a slightly warped accretion disc forms with no more than 3° of 
disc tilt. The stream-disc interaction for the deflected stream forces a large degree of 
material above the orbital plane, increasing the out-of-plane flow drastically. Plotting 
the Ha emissivity in velocity space allows us to compare our results with tomographic 
observations. Deflecting and boosting the stream increases the emissivity in each v z 
slice of the out-of-plane flow by at least three and up to eight orders of magnitude 
compared to the undeflected case. We conclude that a deflected stream is a viable 
mechanism for producing the strong out-of-plane flows seen in the tomographic images 
of U CrB and RS Vul. 
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1 INTRODUCTION 

The recent use of three-dimensional Doppler tomography to 
investigate the Algol binarie s U CrB (|Agafo nov, Sharova fc 



Richards et al. (20101, who noted the similarity in the out- 



Richards, 2009 1 and RS Vul ( Richards, Sharova & Agafonov, 



20101 has unveiled a significant amount of detail regarding 



the accretion process in these systems. As well as confirming 
the presence of a number of features seen in two-dimensional 
tomographic data, the new data have shown an unexpect- 
edly large amount of gas located outside of the orbital plane 
in both systems. U CrB also displays a transient accretion 
disc inclined with respect to the orbital plane, which is not 
expected if mass transfer occurs as a result of gravitational 
and centrifugal forces alone. The companions in short-period 
Algols are rapidly rotating late-type stars that are likely to 
be magnetically active. Magnetic behavior occuring at or 
near the first Lagrange point (LI) may offer an explanation 
for the out-of-plane material. 

The notion that large-scale magnetic effects may dras- 
tically influence the morphology and time-evolution of the 
accretion processes is evidenced by the observations of 



of-plane velocities located at LI and magnetically active re- 
gions on the donor star in RS Vul. Moreover, they found 
evidence of a loop prominence on the cool star close to the 
LI point in the 3D tomogram of RS Vul, which suggests 
that the magnetic field may have deflected the gas stream 
relative to the central plane. Such a deflection would cause 
gas to be introduced into the system with additional (non-z) 
angular momentum that could then propagate toward the 
primary and persist as a tilted disc. This deflection could 
be attributed to a variety of mechanisms. Starspots on the 
donor in the prototype Algol system ft Persei were observed 

~jTMr 



by Olson 



and coronal ejections in eight other Algol 
systems, including RS Vul, have been detected by |White| 
& Marshall (19831. More recently, a large-scale persistent 



coronal loop in Algol was fully resolved by Peterso n et al.| 
(20101, suggesting the possibility of an extended magnetic 
field structure between the two stars. While the likelihood of 
magnetic activity in U CrB and RS Vul is high, the lack of in- 
formation presently available regarding the exact magnetic 
structure of these systems makes it difficult to determine 
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the degree to which magnetic effects are responsible for the 
transport of gas out of the orbital plane. 

Further complicating the mass transfer process in Algols 
is the fact that their short orbital period (< 5 — 6 days) com- 
bined with the close orbit and relatively large radius of the 
primary results in the accretion stream directly impacting 
or grazing the surface of the accreting star. Direct-impact 
systems display a permanent or transient Keplerian disc and 
sometimes a partial disc or localized region characterized by 
an asymmetric mass distribution and sub-Keplerian veloc- 



ities (Richards 1992 Richards & Albright 1999). In these 



systems, the localized region is associated with gas that has 
circled the accreting star and has slowed down because of the 
interaction with the incoming gas stream ( |Richards|[l992l ). 
Gas is expected to leave the orbital plane due to interactions 
at the site where the stream impacts the surface of the pri- 
mary, the stream-disc interaction, and the localized region 
( Richards|2004 l. 

Both U CrB and RS Vul are short-period binaries, and 
both display a number of Ha emission features characteris- 
tic of direct-impact systems. U CrB has been classified as 
an alternating system, since it has been observed in both 



a streamlikc accretion state and a disclike state (Richards 
|fc Albrightj |1999|). The stream state is characterized by 



strong emission along the predicted ballistic trajectory of 
the accretion stream against a lack of circumprimary emis- 
sion. In the disc mode, emission from the accretion disc 
is the most prominent feature. The mechanism responsi- 
ble for the transition between these two states has yet to 
be explained, although hydrodynamic simulations of /3 Per 
I Blondin, Richards & Malinowski, 1995 1 showed that the 



morphology of the accretion flow depended strongly on the 
density of the inflowing gas. The hydrodynamic simulations 



of Richards & Ratliff ( 1998 1 also measured a large rela- 



tive increase in the Ha emissivity of the stream relative to 
the disc when stream density was varied. The similarity of 
the stellar components and orbital parameters of RS Vul to 
those of U CrB (RS Vul = B5V + Gl III-IV, P=4.48 days; 
U CrB = B6V + GO III-IV, P=3.45 days) makes it likely 
that a disc mode is possible in RS Vul, and it has been clas- 



sified as an alternating Algol by Richards & Albright ( 1999 1 



despite having not yet been observed in a disc mode. 

Three-dimensional tomographic imaging of U CrB has 
allowed for a more thorough analysis of the accretion fea- 
tures by enabling the creation of 3D velocity space plots of 
the Ha emission sources. Mapping this data to Cartesian co- 
ordinates is not possible due to the inability to resolve the 
exact Cartesian structure of the system, but the theoreti- 
cal analysis of key features of the system provides reference 
points from which the gas dynamics of the system can be 
inferred. Features commonly used as landmarks include the 
motion of the stellar components, the ballistic trajectory of 
the accretion stream, and the locus of the accretion disc 
as defined by the Keplerian velocities at the surface of the 
primary and its outermost Roche surface. 

The first 3D tomographic analysis of U CrB was per- 
formed by Agafonov, Richards & Sharova (2006) and 



showed the system in a disc mode. An additional series of 
Ha profiles were later analyzed in three dimensions and dis- 
played a stream mode ( [Agafonov et al. | |2009[ ). Notable fea- 
tures of these two modes include: (1) A gas stream transfer- 
ring mass from the donor to the primary star with velocities 



as high as 450 km s _1 in the orbital plane and negligible 
z-velocities. (2) Circumprimary emission associated with an 
accretion disc that displays both positive and negative z- 
velocities (v z = ±420 km s _1 ). The different signs of the z- 
velocities on opposite sides of the primary star suggest that 
the disc may be inclined with respect to the orbital plane. (3) 
A stream-disc interaction region, or "hotspot" (v z — +180 
km s -1 to +540 km s^ 1 ), as well as a star-stream impact 
site (v z — -300 km s _1 to +60 km s _1 ). (4) A localized re- 
gion (v z = +200 km s _1 to +540 km s _1 ). In contrast to 
the two accretion modes displayed by U CrB, RS Vul has 
only been observed in a stream state with a much stronger 
out-of-plane flow. A three-dimensional tomographic analy- 



sis of RS Vul was performed by Richards et al. (20101, and 



displays many of the same features as U CrB. Major differ- 
ences from U CrB include: a truncated gas stream (v z — —85 
km s _1 ) with no indication of a star-stream hotspot, asym- 
metric circumprimary emission with no well-defined disc, a 
localized region diametrically opposite the circumprimary 
emission with opposite z-velocities (suggesting that the gas 
orbiting the accreting star is flowing along a tilted path with 
v z — +15 to —85 km s _1 ), and emission from a magnetically 
active region on the donor that suggests the presence of a 
loop prominence (v z — —150 to 150 km s _1 ). 

The fact that the Ha emission at LI possesses z- 
velocities equal to those of magnetically active regions on 
the donor strongly suggests that the stream is magnetically 
deflected at or near LI. The relatively high velocities asso- 
ciated with the stream also suggest that the stream veloc- 
ity may be boosted to many times the speed that would 
occur via Roche lobe overflow alone. Our goal was to de- 
termine whether a deflected stream is responsible for in- 
troducing a sufficient amount of tilt to the accretion disc, 
and to what extent the deflection influences the presence or 
absence of other out-of-plane accretion features. To investi- 
gate this we performed hydrodynamic simulations of mass 
transfer in RS Vul in which magnetic effects were mimicked 
by manually deflecting and boosting the stream at LI. The 
resulting steady-state solutions were analyzed to determine 
the effect of the deflection on the morphology of the circum- 
stellar material, the presence of out-of-plane flow, and the 
presence or absence of the accretion features noted in the 
tomographic data. The hydrodynamic model we used is de- 
scribed in section 2, results from runs with undeflected and 
deflected streams are given in section 3, an analysis of the 
tilt of the disc with comparisons to analytical estimates is 
given in section 4, and our results summarized in section 5. 



2 COMPUTATIONAL METHODS 
2.1 Grid parameters 

Our simulations were performed with the three-dimensional 
hydrodynamics code VH-1, which was developed by the Nu- 
merical Astrophysics Group at the Virginia Institute for 
Theoretical Astrophysics. VH-1 solves the Euler equations 
for the flow of an inviscid fluid characterized by an adiabatic 
index 7. To solve the system, VH-1 utilizes operator splitting 
to reduce the grid into a number of 1-D segments. The fluid 
equations are solved along each segment using the Piecewise 
Parabolic Method, an explicit high-order Godunov method 
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Figure 1. Yin-yang grid geometry. The figure on the left shows 
the range of one grid, and the figure on the right shows the overlap 
between the two grids. 



tions between the atmosphere and the stream could have a 
non-negligible effect on the structure of the circumstellar gas 
(especially in the case of a rotating and precessing primary), 
but for the purposes of this paper we neglect the structure 
of the atmosphere. Instead, we employ a partially reflecting 
boundary condition that sets the radial velocity of the gas 
at the innermost radial zone to zero every timestep. This 
approximates the gas encountering a shock at the surface 
of the star and accreting shortly afterward. At the outer 
boundary, gas should be allowed to freely exit the grid, so 
we implement a modified zero-gradient boundary condition 
with the added requirement that gas may not be allowed to 
flow inward once it has passed off the grid. 



2.2 Radiative cooling and stream density 



employed on an Eulerian (fixed) grid ( Colella & Woodward 
|1984[ ). The PPM can be used for both 2D and 3D geometries 
without necessitating significant changes to the algorithm 
other than the extension from 2D polar to 3D spherical co- 
ordinates. 

We adapted VH-1 to model our system in spherical co- 
ordinates on a corotating grid located in the reference frame 
of the binary system and centered upon the primary star. 
The inner boundary of the grid coincides with the surface 
of the primary at a radius of Ri, and the outer boundary 
extends to twice the distance from the centre of the primary 
to LI (~ 6.2J?i). Since the system does not vary widely in 
scale between the two stars, we space our zones equally in 
the radial direction. The parameters of U CrB and RS Vul 
are given in Table 1, and were taken from |Agafonov et al.| 
( 20091) andlRichards et al.l (120101). We note that while we 



use the parameters corresponding to RS Vul in our simula- 
tions, both the orbital parameters and stellar components 
of U CrB are similar to those of RS Vul. This enables us to 
interpret our results in the context of either system. 

To avoid grid convergence and the coordinate singulari- 
ties intrinsic to spherical coordinates, we employ a yin-yang 



overset grid (Kageyama & Sato 20041. The yin-yang grid 



consists of two separate spherical polar grids oriented 90° 
with respect to each other, so that the polar region of one 
grid aligns with the equatorial region of the other. Each 
grid extends ±7r/4 azimuthally, and ±37r/4 in the polar di- 
rection (Figure [TJ . This grid configuration allows for rea- 
sonable limits on the Courant number and ensures that our 
runtimes do not become excessive. As a result, we are able 
to extend our outer boundary to a sufficiently large radius 
without imposing limits on the the azimuthal extent of the 
grid. We define our coordinates such that the positive x-axis 
of the yin grid lies along the binary axis and points toward 
the secondary. The positive y-axis points in the direction of 
motion of the secondary, and the z-axis is perpendicular to 
the orbital plane. Each piece of the grid has dimensions of 
288 x 144 x 432, which results in an angular resolution of 
A0 ~ 0.625° and Ar/Rx = 0.018 at the inner boundary. 

Since the inner boundary corresponds to the surface 
of the accreting star, we must choose a boundary condition 
that allows gas to accrete while still adequately modeling the 
effects of the stellar atmosphere. The fact that the stream 
penetrates only a small portion of the stellar atmosphere 
( Richards|l992| ~ 10" 6 7?i for Algol, ) allows us to approxi- 
mate the surface of the primary as a solid surface. Interac- 



The 2D simulations of Algol by Blondin et al. (1995) found 



that the morphology of the circumprimary flow depended 
strongly on the initial stream density due to the effects of 
radiative cooling. The stream density's influence on the flow 
is seen primarily in the stream-disc interaction. A denser 
stream will cool more rapidly at the shock, negating the 
effects of heating within the shock. This results in a flow 
that behaves similarly to a compressible isothermal gas. A 
lower-density stream will not cool to the same extent, and 
will behave adiabatically. 

To investigate the relative effect of cooling in U CrB 
and RS Vul, we compare the cooling time at the shock to a 
characteristic dynamical time. The dynamical time can be 
found using the gas speed at the outer edge of the disc where 
the stream-disc shock is located. Assuming a Keplerian ve- 
locity at the outer edge of the disc (~ 350 km s~l) allows 
us to calculate a dynamical time of: 



GMi 



u Kep 



(1) 



In both U CrB and RS Vul, the dynamic timescale is 
calculated to be 0.04P or f,. 

If a cooling function is known, we can calculate the char- 
acteristic cooling time at the shock to be 



t c = 



kT 



n s A(T) 



(2) 



Here, n s is the electron number density in the post- 
shock flow, and A is the cooling function evaluated at the 
temperature of the shocked gas, T. We calculate A(T) with 
the cooling function for a low-density plasma as given by 
|Cox &: Daltabuit| ( |1971| ). From the Rankine-Hugoniot jump 
conditions, we can calculate the post-shock temperature to 
be T = 3mn^/16fc, where m is the average mass per parti- 
cle. We assume a gas comprised largely of ionized Hydrogen, 
in which case rh ~ mjj/2. We assume the velocity of the 
shocked gas, v s h to be on the order of the Keplerian velocity 
at the shock site. For both U CrB and RS Vul, we calculate 
T s h to be approximately 1 x 10 6 K, which results in a cooling 
coefficient of A = 8 x 10~ 23 ergs cm 3 s _1 . An analytical es- 
timate of the density at the shock site is difficult to obtain, 
so we note that our simulations show densities in the shock 
region on the order ofp = 0.1pLi - Converting this density to 
physical units requires a measurement of the mass transfer 
rate. Since data on the exact mass transfer rates for U CrB 
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Table 1. Parameters for U CrB and RS Vul 



Parameter 


U CrB 


RS Vul 


Mass of Primary, Mi 


4.8OM 


6.59M Q 


Mass of Secondary, M2 


1.40M© 


1.76M 


Effective Temperature of Primary, Ti 


14000 K 


16500 K 


Effective Temperature of Secondary, T2 


5250 K 


6000 K 


Radius of Primary, i?i 


3.OR 


4.7^ 


Radius of Secondary, R2 


4.6R Q 


5.8R 


Binary Separation, a 


1.23 X 10 12 cm 


1.61 x 10 12 cm 


Orbital Period, P or t, 


3.452 days 


4.478 days 



and RS Vul are not available, we use the l ower limit of th e 
transfer rate in Algol, M w 10 _11 M Q yr" 1 |Richards|l992 1. 



Comparing the characteristic dynamical time to the 
cooling time, we find that ta/t c ~ 2.4. This is a conser- 
vative estimate, as we have used a small value of the mass 
transfer rate to calculate the cooling time. Any increase in 
the mass transfer rate will serve to decrease the cooling time 
further, since t c varies as 1/M. In general, the mass transfer 
rates in short-period Algols vary widely between systems, 
and are not likely to remain constant in ti me. For Algol, M 



could range from 10 _li to 1O" M M yr" 1 (Richards 19921 



With these factors taken into account, it is likely that cool- 
ing at the shock is sufficiently rapid to counter the heating 
that occurs in the stream-disc interaction. This rapid cooling 
ensures that the shock-heated gas will quickly fall to some 
minimum temperature below 10 4 K, at which point the re- 
combination of H will diminish the cooling rate. The cooled 
gas will be heated by the strong photoionizing UV flux from 
the primary, which is sufficient to maintain the temperature 
at a value close to 10 4 K. Taking this information into ac- 
count, we assume a constant temperature of 11, 000K in our 
simulations. This value falls midway between the effective 
temperatures of the primary and secondary in both U CrB 
and RS Vul. This allows us to expedite our calculations by 
employing an isothermal equation of state (7 = 1). 

2.3 Initial conditions at LI 



Lubow & Shu ( 1975 hereafter LS) conducted an extensive 



analysis of the gas dynamics at the LI point, and we use 
their prescription to determine the initial conditions of the 
accretion stream. LS parametrized the flow in terms of the 
ratio of isothermal sound speed to the system's orbital speed: 
£ = c a /bja = (fcT/m) 1//2 , for which we obtain e = 0.035. This 
gives a local sound speed of 9.5 km s _1 at LI. LS calculate 
that the stream will be accelerated to supersonic speeds as 
it travels through LI, so we assume an unboosted speed of 
Mach 1. The angle the stream makes with the binary axis 
can also be calculated using e, and for RS Vul we obtain 
6 = 20.4°. Assuming an ideal gas, the pressure at LI should 
be initialized to a value on the order of pe 2 . 

For this work we require control over the direction and 
magnitude of the gas velocity at LI in order to manually de- 
flect the stream. To achieve this, we calculate the position 
of LI on the grid and hold the stream parameters in the 
corresponding zones fixed at the values given above. This 
eliminates the need for the donor star to be modeled hydro- 
dynamically, and we remove it from the grid while retaining 
its term in the gravitational potential. A small amount of 



flow into the region corresponding to the secondary occurs, 
but it remains confined to the volume that the secondary 
would normally occupy and does not affect the behavior of 
the accretion flow. In addition to giving us more control over 
the stream at LI, this strategy has the benefit of allowing 
for a more economic grid size, since the full secondary need 
not be present on the grid. 

LS calculate the width of the stream at LI to be on 
the order of ea, which requires an angular resolution of at 
least A(j> = 0.022. The angular resolution of our grid is small 
enough to accomodate this with two zones; each has an an- 
gular extent of A<f> — 0.011. Symmetry requires the gas 
stream at LI to extend equal distances above and below 
the orbital plane as well as an equal distance to either side 
of the binary axis. This necessitates the use of four zones to 
represent LI. In reality the density profile of the gas at LI is 
Gaussian, but we lack the resolution to resolve this. Instead 
we hold the density of the gas fixed to a constant at our LI 
zones. Ultimately, our system exhibits behavior similar to 



that of Blondin et al. ( 1995 I in that the morphology of the 



accretion flow is relatively insensitive to the geometry of the 
zones at LI. 

Our simulations do not assume any gas to be initially 
present around the primary, and the system requires many 
orbital periods to allow the full formation of a disc. A steady- 
state solution in which a disc is present and has stopped 
gaining mass (indicated by M acc = Mli) requires 40 — 50 
Port according to preliminary low-resolution runs. To expe- 
dite the disc formation, we multiply the density in the disc 
by a constant factor after two orbital periods have elapsed. 
Since the mass accretion rate is proportional to the density 
at the inner boundary, this results in the mass accretion rate 
being boosted to a value closer to that at which it stabilizes. 



3 RESULTS 

3.1 Undeflected stream 

To ensure that our model behaved appropriately in the ab- 
sence of a deflected stream, we performed a fiducial run with 
the stream left unmodified (0 = 0°, U stream 

= 1M). Figure 

[2] shows the stream following the ballistic trajectory as it 
falls toward the primary and narrowing as it approaches the 
surface. The stream impacts the surface after approximately 
O.lPorb- We allow gas to accrete until 2P or t, during which 
a thin accretion disc begins to form. At 2P or b, we boost 
the disc density. The boost instigates some radial density 
oscillations due to increased pressure along the density dis- 
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Figure 3. Density vs. radius, averaged over 4> within the orbital 
plane at 18P or t- Densities are normalized to 1 at LI. The solid line 
depicts the density for the undeflected stream, while the dashed 
line corresponds to a stream tilted to 45° and boosted to Mach 
30. LI is located at 0.63a, and the surface of the primary at 0.20a. 
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Figure 4. Tilt of the disc vs. radius, calculated at 18P or b by 
summing the angular momentum within the disc over cf> for each 
value of r. The solid line depicts tilt for the undeflected stream, 
while the dashed line corresponds to a stream tilted to 45° and 
boosted to Mach 30. 



continuity that occurs at the radius of the boosted region. 
These oscillations propagate with the flow around the disc, 
but fade within the next few orbital periods. Following the 
density boost, the radius of the disc continues to increase 
for several orbital periods, and a steady state is achieved at 
roughly 15P orb . 

When the system has stabilized, we measure the disc 
to have an outer radius of 0.34a (where a is the binary sep- 
aration as measured from the center of the accreting star). 
Examining the density profile averaged over 4> (Figure |3| 
shows that the density of the innermost portion of the disc 
(~ 0.22a) is of the same order of magnitude as the den- 
sity at LI (~ 0.63a). The density decreases sharply to just 
above ~ 10 _4 pLi at the outer edge of the disc. The disc is 
symmetric about the orbital plane to approximately ±0.03a 
above and below the plane. Beyond this region the density 
decreases rapidly, and is measured to be less than 10 _4 pLi- 

The velocities within the disc are oriented largely within 
the plane. As the gas accelerates from LI, its speed increases 
from ~ 10 km s _1 to a final velocity of ~ 240 km s _1 be- 
fore encountering the stream-disc shock. (All velocities are 
measured with respect to a corotating frame.) After cross- 
ing the shock and entering the higher-density regions of the 
disc, the gas assumes velocities which fall within the range 
of 240 — 440 km s , with small velocities corresponding to 
the outer edge of the disc, and higher velocities located near 
the surface of the accretor. The velocities within the disc re- 
main relatively constant in direction and are sub-Keplerian 
throughout (60 — 80 percent at the outer and inner edges 
of the disc, respectively). These outer velocities are less uni- 
form near the localized region and stream-disc interaction 
site, and reach values up to 300 km s _1 in the vicinity of 
the stream-disc shock. 

To quantify the angle with which the disc is inclined, 
we integrate the angular momentum present in the disc and 
measure its angle with respect to the z-axis. We limit our 
integration to the region occupied by the high-density por- 
tions of the disc, using Figure [3] as a guideline. The result is 
summed angular momentum values for each value of (r, (j>) 
in the orbital plane. Each value of r defines a ring within 



the disc for which we are able to calculate the tilt. The sign 
of the tilt is set equal to the sign of the x-component of 
the net angular momentum. This has the effect of assigning 
positive tilt angles to rings that are tilted toward the donor, 
and negative angles for rings that are tilted away. For the 
undeflected stream, the disc is largely symmetric about the 
orbital plane and remains essentially untilted over time. Ex- 
amining the tilt as a function of radius at 18P or b (Figure 
Q shows that the high-density portions of the disc from the 
inner boundary to ~ 0.3a only possess tilt on the order of 
0.1°. This is a fraction of the angular resolution of by one 
grid zone. The tilt varies gradually as the radius increases, 
and the sign changes three times along the radius, indicat- 
ing that some warping occurs about the orbital plane. The 
tilt becomes slightly larger in the low-density portions of 
the disc and reaches a maximum value of 0.5°. This tilt 
is smaller than the amount able to be resolved with one 
zone. The presence of this tilt is likely a result of a compar- 
atively large amount of non-z angular momentum created 
in the stream-disc interaction through the splashing of low- 
density gas. Turbulence in the gas flowing along the outer 
low-density edge of the disc may also contribute to this mea- 
surement. 

The disc's orientation with respect to the orbital plane 
can be found by examining a plot of the x and y components 
of the angular momentum. In cartesian space, the projec- 
tion of the angular momentum vector into the orbital plane 
points toward the direction that the disc is azimuthally in- 
clined. (This is analogous to treating the disc as a rigid body 
with a rod through its axis of rotation. The projection of the 
angular momentum vector into the orbital plane is the direc- 
tion you would push the rod to tilt the disc.) Figure [5] shows 
a J y vs. J x plot for both runs, with J x and J y components 
shown as a fraction of the total angular momentum. The 
data points for the undeflected stream lie clustered around 
the origin, with magnitudes of less than 1 percent of the total 
angular momentum. This supports the qualitative observa- 
tion that the disc lies almost entirely in the orbital plane 
throughout the length of the simulation. The disc appears 
to precess at a period of P or b, which is a by-product of mea- 
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Figure 2. Top and side views of the primary star and and disc for an undefiected stream after 
18P or f,. The density scale is normalized to 1 at LI. The top view depicts the density on a slice 
through the orbital plane. Density isosurfaces in the side view are taken at p = 1 X 10 — 1 ,5 X 
10 -2 , 1 X 10 -2 , 5 X 10 — 3 , and 1 X 10 — 3 with the transparency increased by a factor of 0.5 for each 
isosurface. 
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Figure 5. Components of the angular momentum of the disc for 
the undefiected stream (circles) and the deflected and boosted 
stream (crosses). Data points were taken at intervals of ~ 0.2P or ;,. 
The undefiected stream results in an disc that remains largely in 
the orbital plane, with minor variations over time that remain 
symmetric about the origin. The deflected and boosted stream 
produces a disc biased toward +J X due to the stream impacting 
the far side of the disc from below and tilting it toward the donor 
star. 

surements of the angular momentum taken from within a 
co-rotating reference frame. 

3.2 Deflected stream 

To simulate a system with a deflected and a boosted veloc- 
ity, we modified the stream parameters at the zones corre- 
sponding to LI by directing the stream velocity 45° below 
the orbital plane and boosting it to thirty times the local 
sound speed; this is comparable to gas flows within magnetic 
loops. These values were chosen to maximize the likelihood 
that a tilted disc would form. In the early stages of evolution 
(Figure [fjjl, the stream falls along the predicted ballistic tra- 
jectory, which extends below the orbital plane to just below 
the bottom of the accreting star. As it nears the primary, 
it is pulled back toward the orbital plane and wraps around 
the accretor. The boosted velocity is large enough to cause 
the stream to miss the surface of the star, resulting in the 
absence of a star-stream collision and a gas stream that ex- 
tends above the orbital plane to nearly 2a. The stream then 



circles the accretor and interferes with the gas stream from 
the secondary as an accretion disc begins to form near the 
orbital plane. The flow during the first few orbital periods 
displays a large degree of out-of-plane flow as the system 
approaches a steady state. At 2P or (, the density within the 
disc is boosted, inducing some minor radial oscillations as 
in the undefiected case. As the disc gains mass, the stream- 
disc interaction continues to drive gas upward above the 
orbital plane. The system stabilizes over the next 15P or b, 
after which the accretion disc remains present along with a 
large plume of gas deflected above the orbital plane by the 
stream-disc interaction. 

Investigating the morphology of the accretion disc 
shows some radial enlargement and warping due to the de- 
flected stream. Averaging the density in the central plane 
over cf> (Figure[3]| shows that the density falls off more slowly 
than in the undefiected case, reaching 10 _3 pLi at 0.5a. The 
inner density of the disc remains on the same order of mag- 
nitude as in the undefiected case, and is slightly larger in 
radius. The top view of Figure [6] at 18P or (, shows that the 
disc is also asymmetric, with a slightly enlarged portion in 
the upper right portion of the disc. This large density region 
originates as the disc forms and propagates around the disc 
with time. It does not appear to originate due to the artifi- 
cial density boost that occurs at 2P or (,. The vertical extent 
of the disc is similar to that of the undefiected case, neglect- 
ing slight enlargement above and below the orbital plane 
due to warping that occurs in the outer edges of the disc. 
The majority of the disc still lies in the orbital plane, but 
with lower density gas (~ O.lpLi) deflected above the plane 
due to the stream-disc interaction. This gas is pulled back 
into the plane, where further dissipation of angular momen- 
tum due to interactions with the disc causes a portion of the 
material to join the disc and eventually accrete. 

The deflected stream begins at a speed of 285 km s~ 
(Mach 30). As the circumprimary material begins to build 
up into a disc, gas is blown upward in the collision with z- 
velocities of up to 210 km s _1 . Within the disc we observe 
a circular orbit with velocities that remain sub-Keplerian in 
the same range as the undefiected case (240 — 440 km s" 1 
for the high-density portion, 150 — 250 km s _1 for the low- 
density portion). We also measure z- velocities on the order 
of ± 10 km s _1 due to the warping and tilting in the disc. 

Examining the tilt within the disc as a function of radius 
after 18P or b has elapsed (Figure [fl) shows that the warping 
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Figure 6. Top and side views of the primary star and disc for a stream boosted to Mach 30 
and deflected 45° below the orbital plane. The density scale is normalized to 1 at LI. Density 
isosurfaces are taken at p = 1 X 10" 1 , 5 X 10~ 2 , 1 X 10~ 2 , 5 X 10 -3 , 1 X 10 -3 with the transparency 
increased by a factor of 0.5 for each isosurface. Due to asymmetry about the orbital plane, density 
isosurfaces are depicted in the top view, as opposed to a slice along the plane. 



within the disc is more pronounced than in the undefiected 
stream. The inner portion of the disc is tilted at ~ 1.5° 
near the the surface of the primary, and increases to ~ 3° 
at 0.35a. Warping is indicated by the variations in tilt angle 
with increasing radius. The tilt continues to increase at radii 
greater than 0.35a, mostly due to the high z- velocities in the 



material blown upward in the stream-disc interaction. The 
tilt angle calculated for this low-density gas reaches values 
as high as 50° at the outer edge of the disc. All the tilt mea- 
surements for this timestep result in positive values, indi- 
cating that the disc is tilted consistently in the +x direction 
(toward the secondary). 
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The orientation of the disc changes noticeably over time 
when we deflect and boost the stream. The plot of J y vs. J x 
(Figure [5| shows a much wider spread of values than that of 
the undeflected stream. The fraction of angular momentum 
that is distributed in the x and y components remains be- 
tween 2 and 3 percent at all times, and grows to 5 percent 
when the disc is oriented toward the donor star. When this 
occurs, the effects of the stream's impact on the disc are en- 
hanced due to the upward force of the stream on the already 
upwardly-moving gas in the disc. As in the undeflected case, 
the apparent precession is due to measurements from within 
the corotating reference frame. 

3.3 Emissivity 

Additional information regarding the out-of-plane velocities 
and accretion features are more easily visible if we mimic the 
tomographic results with a pseudo-emissivity plot in veloc- 
ity space. This w as first done for Algol systems by Richards 
& Ratliff (1998). In our simulations, emissivity scales as p , 



and lacks any dependence on the temperature or ionization 
state of the gas. We slice our velocity space plots into slabs 
along v z , sum the emissivity along v z within each slab, and 
scale the emissivity to the maximum value observed in the 
slab. This maximum value is also listed in terms of the max- 
imum emissivity on the entire grid. The velocities measured 
in our hydrodynamic simulations take place in a corotat- 
ing frame, and we transform these into an inertial frame for 
the sake of comparison with tomographic data. Since we use 
the inner boundary to represent the primary and we have 
removed the secondary from the grid, emission features asso- 
ciated with stellar activity occuring on the surface of either 
star (e.g. circumprimary emission, chromospheric emission 
from the secondary) are absent from the plot. As a result, 
the plots appear more sparse than the original tomographic 



data presented in Agafonov et al. ( 2009 1 and Richards et 
|al.| ( |2010[ ) . We also note that the strong density dependence 
causes the emissivity in the central plane to appear many 
orders of magnitude greater than that of other slices. 

The emissivity plots for the undeflected stream (Fig- 
ure [7| clearly show the presence of the disc in the central 



slice (v z 



-30 to 30 km s ). Since the inner portion of 



the disc circles the primary with a higher velocity than the 
outer portion, the disc is turned inside-out on the velocity 
space plot. The bulk of the disc possesses velocities that fall 
within the Keplerian predictions (400-515 km s _1 ), although 
we also measure sub-Keplerian velocities as small as 370 km 
s _1 . The accretion stream is also visible in the central slice, 
beginning near v x = km s -1 and extending further in the 
-x direction as it gains speed and falls along the binary axis. 
The slightly stronger emissivity region located along the in- 
side region of the "ring" at roughly (v x ,v y ) — (—250,250) 
km s" 1 corresponds to the stream-disc interaction. Emission 
from this feature extends up to v z = 30 km s _1 . We also note 
that the emissivity is slightly enhanced near (±350, —350) 
km s _1 due to the overlap between the yin and yang grids. 

Some emission is present in the slab corresponding to 
v z = 30 to 90 km s _1 , and is located primarily within the lo- 
cus of the accretion disc. The intensities of this emission are 
seven orders of magnitude below the maximum intensity on 
the grid, and distinct emission features are not immediately 
apparent in this slice. For gas wtih downward velocities of 



v z = —30 to —90 km s _1 , some emission is present in the 
locus of the disc again, though at intensities ten orders of 
magnitude below the maximum grid intensity. For veloci- 
ties greater than v z = ±90 km s _1 , the emission lacks any 
distinct features and occurs at intensities on the order of 
10 -13 — 10~ 14 of the maximum intensity. 

The emissivity plots of the tilted disc (Figure |8| show 
much higher emissivity measurements in slices above and 
below v z — 0. The disc is again the dominant feature in the 
central slice, and is larger in radius than in the undeflected 
case, as expected from the density profiles. This results in an 
emissivity region that extends inward to a slightly greater 
extent. The majority of the velocities lie within the Keple- 
rian predictions (375—515 km s" 1 ), with some gas possessing 
velocities as low as 350 km s" 1 . Since the stream has been 
deflected below the orbital plane, it appears in the upper 
left quadrant of the — v z slices. The stream-disc interaction 
is visible in the slices corresponding to positive z-velocities, 
and appears most prominently near (v x ,v y ) — (—450,150) 
km s _1 in the slice of v z = 30 to 90 km s" 1 . The high v z 
material on the far side of the disc is also apparent in the 
emissivity plots. As it rises and wraps around the secondary 
(with +v z ), it can be seen in the lower left quadrant. As it 
falls back towards the orbital plane (— v z ), it shifts toward 
the lower right quadrant, indicating motion back toward LI. 
The gas collides with the orbital plane again, and a portion 
is driven below the orbital plane once again. This is visible 
in the +v z slices as a feature in the upper right quadrant. 

As a whole, the emissivity for the deflected and boosted 
stream reveals values above and below v z = with intensi- 
ties much closer to the values measured in the central slice. 
This indicates a much stronger out-of-plane flow produced 
by the deflected and boosted stream. It should also be noted 
that the maximum intensity on the grid is nearly five times 
larger for the deflected and boosted stream than for the un- 
deflected stream. This is to be expected, since the densities 
within the disc are slightly larger when the stream was de- 
flected, and any changes in density will appear in the emis- 
sivity as p 2 . Examining the emissivity in the early stages 
of our run (2P OT i, just before the density is boosted) shows 
the intensity more uniformly distributed throughout the v z 
slices. No notable differences in the structure of the emissiv- 
ity features are seen. 



4 ANALYSIS 

The main features visible on our emissivity plots for both 
the undeflected and deflected streams are the disc emission, 
accretion stream, and stream-disc interaction. In the unde- 
flected case we also measure non-localized "clouds" with z- 
velocities up to ±210 km s" 1 , but with magnitudes smaller 
than the emission features located in the central plane by 
a minimum of eight orders of magnitude. For the deflected 
stream, the main sources of emissivity are the disc and the 
out-of-plane gas located on the far side of the secondary 
(extending from v z = +210 km s _1 to v z = —150 km s^ 1 ). 
The disc structure is similar to that of the undeflected case, 
but displays some radial enlargement and slight warping. 
The out-of-plane gas appears with much stronger intensi- 
ties relative to the gas in the orbital plane, and measures 



10" 



10" 



as opposed to 10 — 10 
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Figure 7. Emissivity for RS Vul with an undeflected stream after 18P or (,. The emission intensity 
on each plot is normalized to the maximum within the individual slice. The value of each maximum 
is also listed in terms of the maximum intensity on the entire grid. 
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Figure 8. Emissivity for RS Vul with an accretion stream deflected below the orbital plane by 45° 
and boosted to Mach 30 at 18P oro . As in the previous figure, the emission intensity is normalized 
to the maximum value within each individual slice, and the value of each maximum is listed in 
terms of the maximum intensity on the entire grid. 
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The early stages of our deflected stream simulation are 
comparable to the tomographic data of RS Vul ( |Richards| 
et al.pOlO l. The emission seen in the tomographic data of 
RS Vul is not confined to the Keplerian locus and is found 
primarily along the stream, near the mass donor and ac- 
cretor, and in several other localized high-velocity regions. 
In our simulations, prior to the formation of the accretion 
disc (i.e. before the density boost occurs at ~ 2P or b) we ob- 
serve a similar structure. The stream is visible in slices with 
negative velocities as low as v z — —210 km s" 1 , and the 
high- velocity features corresponding to the stream material 
wrapping around the primary are visible from v z — —150 
km s _1 to v z = +210 km s _1 . As in the tomographic data, 
the star-stream interaction is not visible, since the deflected 
stream misses the star entirely. Emission from the Keple- 
rian locus is non-uniform and occurs at intensities lower than 
those corresponding to the stream by an order of magnitude. 

Despite similarities with the morphology of our simu- 
lations, the tomographic data from RS Vul show a more 
feature-rich display, with out-of-plane velocities as high as 
v z = ±480 km s _1 . These are much larger than the ±210 
km s _1 maximum we measure in our emissivity plots, but 
the magnitude of the tomogram velocities may be some- 
what exaggerated due to the data's relatively low resolution 
in the z-direction. Another important difference from our 
model is that the intensities of the emission features in the 
tomograms remain on the order of those in the central plane 
even at very high values of v z . The primary explanation for 
our lower intensity measurements is that the emission of 
low-density regions is suppressed by our use of an isother- 
mal gas. Since the emissivity calculation depends strictly on 
density in the isothermal case, the magnitude of the emis- 
sion in the simulation data appears smaller than it other- 
wise would in regions containing hot low-density gas. Since 
we have neglected the influence on the ionization state and 
temperature of the gas in our emissivity plots, we expect our 
intensities to fall drastically in the lower-density portion of 
the flow. In the early stages of the run, this results in the 
maximum intensities occurring at the location of the stream. 



The intensity elsewhere is smaller by a factor of 10 



-10" 



The late stages of our simulation display the accretion 
disc in a steady-state with a morphology more akin to the 



tomographic data of U CrB in its disc mode (Agafonov et 
aL"||2009| . This tomographic data shows more concentrated 



emission in the Keplerian locus than the tomograms of RS 
Vul, as well as both a star-stream and stream-disc interac- 
tion region. Our late-stage emissivity plots are also domi- 
nated by emission from the Keplerian locus due to the pres- 
ence of a disc, but the material driven away from the orbital 
plane in the stream-disc interaction remains visible. Due to 
the strong density dependence of the emissivity, the disc 
dominates the emissivity plots. Intensities elsewhere range 



from 10 



10 below the maximum measurement on the 



central slice. 

Our measurements of the tilt angle of the disc show that 
deflecting and boosting the stream introduces warping of up 
to 3°, but no significant tilt with respect to the orbital plane. 
Alternately, the lower-density gas blown upward above the 
orbital plane does display significantly larger tilt that mea- 
sures up to 50° at the outer edge of the disc. Obtaining a 
tilt measurement from the tomographic data for comparison 
is difficult, particularly for RS Vul, which does not display 



a well-defined accretion disc. In U CrB, z-velocities of ±420 
km s" 1 were observed within the predicted locus of the Ke- 
plerian accretion disc. A tilted disc was inferred from the 
observation that emission in tomogram slices corresponding 
to negative z-velocities was stronger on the left side of the 
Vy — v x velocity-space plots, while positive z-velocity slices 
showed stronger emission in the lower-right side of the v v —v x 
plot. Quantifying this tilt by assuming a z-velocity of 420 km 
s _1 and Keplerian velocities at the surface of the primary 
gives a tilt angle of 37.2°. Performing the same calculation 
at the outer Keplerian orbit produces a tilt of 56.1°, which 
is comparable to the tilt measured in our simulations for the 
plume of gas on the far side of the accretor. It should also 
be noted that the v z values are uncorrected for the lower 



resolution in the v z direction. Agafonov et al. (2006) claim 



the correction factor will be small, but any correction factor 
will serve to reduce the predicted tilt. Assuming a factor 
as large as 0.5 would reduce the maximum z-velocities to 
v z = ±240 km s _ , which is close to the maximum v z we 
measure. This correction results in tilts of 20.8° and 36.6° 
at the inner and outer radii. Our simulations indicate that 
it is unlikely that the disc will be uniformly tilted, since 
the stream-disc interaction acts primarily on the outer edge 
of the disc. The tomographic data that depicts asymmetric 
intensities measured on opposite sides of the accreting star 
could be interpreted as warping in the disc, as opposed to 
a uniformly inclined disc. A large part of the out-of-plane 
flow in our simulations also occurs within the Keplerian lo- 
cus, and could contribute to this measurement. 

The primary difference between the deflected and undc- 
flected stream in our simulations is that the modified stream 
enables large amounts of flow above and below the orbital 
plane. Boosting the stream allows gas to circle the primary 
without striking its surface, which in turn means the sys- 
tem retains angular momentum that would normally be de- 
stroyed in the stream-disc and star-stream collisions. The 
additional (non-z) angular momentum manifests in the out- 
of-plane flow, and is eventually dissipated in the accretion 
disc. The disc's dimensions appear larger and the densities 
slightly higher due to the absence of the star-stream interac- 
tion as an angular momentum dissipation mechanism. The 
bulk of the disc remains primarily within the orbital plane, 
although it is more dense and more uniform than the disc 
shown in the tomographic measurements of U CrB. We can 
reconcile the differences between the simulated and observed 
discs by noting that our isothermal approximation results in 
a gas that is more easily compressible. This will lead to a 
more compact accretion disc and the tendency for less gas 
to be blown away in the interaction regions. It is also possi- 
bile that the tomographic images of U CrB in a disc mode 
show the early stages of the disc, in which case the relative 
intensity of the disc would be lower, and the disc would have 
a more nonuniform structure. 

We see no evidence of a transition between a stream 
and disc accretion mode in our simulations. This is likely a 
result of our use of a constant mass transfer rate through 
LI, and potentially a byproduct of our isothermal equation 
of state. Examining our emissivity plots at 2P or t shows a 
morphology very similar to that at 18P or &, but with the 
central slice intensity 1 — 2 orders of magnitude higher than 
that on the remainder of the grid. While the disc is present, 
it is not fully formed, and appears less uniform than late 
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in the run. A reduced rate of mass transfer could offer the 
disc material a chance to accrete without being replenished 
by material from the stream. A non-isothermal gas would 
contribute to this effect by resulting in a less compressible 
flow that would impede the formation of a stable disc. Dur- 
ing this time, we would expect to observe the system in its 
stream phase. Since the time necessary for a full disc to 
form is large (preliminary runs showed Mli = M acc was 
achieved after 40 — 50P or t), variations in the mass trans- 
fer rate occurring across smaller timescales would prevent 
the full formation of a stable disc. Our use of an isothermal 
gas also increases the compressibility, and subsequently the 
density of our disc. Relaxing the isothermal approximation 
could result in a lower-density, less uniform disc. 



5 CONCLUSIONS 

Our simulations of RS Vul have shown that boosting the ve- 
locity of the accretion stream and deflecting it below the or- 
bital plane at Lf introduces significant changes in the struc- 
ture of the accretion flow. Using an isothermal equation of 
state, we simulate the formation of an accretion disc with 
both an undeflected stream and a stream boosted to Mach 
30 and deflected by 45°. In the undeflected case, we mea- 
sure the formation of an accretion disc that is symmetric 
about the orbital plane and displays less than one angular 
grid zone (0.5°) of warping. Modifying the stream causes 
some radial enlargement of the disc, along with warping as 
large as 3°. In addition, a low-density plume of gas is ejected 
above the orbital plane on the far side of the accreting star 
with densities on the order of I0~ 2 pLi — 10 _3 pLi and tilt 
angles as high as 50°. 

The intensities we measure in our pseudo-emissivity 
plots are strongest in the orbital plane due to our use of 
an isothermal gas, which causes the emissivity to depend 
solely on p . The increased out of plane flow that occurs 
due to the deflected stream increases the emissivity above 
and below the central plane by a minimum of five orders 
of magnitude. We anticipate that relaxing the isothermal 
condition would cause higher emissivities to be measured 
outside of the central plane, as in the original tomographic 
data. 

Our results are consistent with the hypothesis that 
stream deflection is responsible for the "tilted" out-of-plane 
material that is measured in the tomograms of U CrB, and 
also for the extreme out-of-plane flow observed in RS Vul. 
We have shown that a strong velocity boost and deflection 
at Lf are sufficient to change the direct-impact nature of 
the system, which increases the amount of out-of-plane flow 
drastically. We expect that a variable mass transfer rate and 
potentially more extreme boosting would only serve to en- 
hance the features observed in the out-of-plane flow. 



award. The author would like to thank Dr. John Blondin for 
his assistance with this work, and the anonymous referee for 
the helpful suggestions regarding the manuscript. 



REFERENCES 

Agafonov, M. I., Richards, M.T., Sharova, O. I. 2006, ApJ, 
625, 1547 

Agafonov, M. I., Sharova, O. I., Richards, M.T. 2009, ApJ, 
690, 1730 

Blondin, J. M., Richards, M. T., Malinowski, M. L. 1995, 
ApJ, 445, 939 

Colella, P., Woodward, P. R. 1984, J. Comput. Phys., 54, 
174 

Cox, D. P., Daltabuit, E. 1971, ApJ, 167, 113 
Kageyama, A., Sato, T. 2004, GGG, 5, 9005 
Lubow, S. H., Shu F. H. 1975, ApJ, 198, 383 
Olson, E. C. 1981, ApJ, 250, 704 

Peterson, W. M., Mutel, R. L., Giidel, M., et al. 2010, Nat, 
463, 207 

Richards, M.T. 1992, ApJ, 387, 329 
Richards, M. T., Albright, G. E. 1999, ApJS, 123, 537 
Richards, M. T., Ratliff M. A. 1998, ApJ, 493, 326 
Richards, M. T., 2004, AN, 325, 229 

Richards, M. T., Sharova, O. I., Agafonov, M. I. 20f0, ApJ 
720, 996 

White, N. E., Marshall, F. E. f983, ApJ 268, f 17 



6 ACKNOWLEDGMENTS 

This research was supported by an allocation of advanced 
computing resources provided by the National Science Foun- 
dation. The computations were performed on Kraken at the 
National Institute for Computational Sciences. This research 
was also supported in part by a North Carolina Space Grant 



© 0000 RAS, MNRAS 000, 000-000 



